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Abstract 

We study the nucleon-pion-state contribution to the QCD two-point function of standard nucleon 
interpolating fields. For sufficiently small quark masses these two-particle states are expected to 
have a smaller total energy than the single-particle excited states. We calculate the nucleon- 
pion-state contribution to leading order in chiral perturbation theory. Both parity channels are 
considered. We hnd the nucleon-pion-state contribution to be small, contributing at the few percent 
level to the effective mass in the positive parity channel. 
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I. INTRODUCTION 


Lattice QCD has made enormous progress over the last years due to computational and 
algorithmic advances [T] . This has led to signihcantly improved lattice calculations of many 
low-energy observables. Present-day unquenched lattice calculations are performed with 
quark masses close to or at their physical value [2H5]. Uncertainties associated with the 
chiral extrapolation are essentially eliminated in these simulations. 

However, dynamical lattice QCD with small quark masses may face new problems. One 
feature of unquenched lattice simulations is the presence of multi-particle states in the 
correlation functions measured to obtain observables. With the up and down quark masses 
getting closer to their physical values one expects multi-particle states with additional pions 
to become a significant excited-state contamination in many correlation functions. As a 
simple example consider the two-point function C{t) of a nucleon interpolating field, as it is 
measured to extract the nucleon mass From the spectral decomposition the two-point 
function in a hnite spatial volume, projected to zero momentum, is a sum of exponentials, 

C'(t) = -t-... . (1.1) 

The first exponential provides the exponential decay with the nucleon mass, Mq = M^y. All 
the other terms stem from states with the same quantum numbers as the nucleon, either 
genuine single-particle excited states or multi-particle states. For sufficiently small pion 
masses one expects a nucleon-pion state to be the state with lowest total energy next to the 
ground state. Mi For symmetry reasons the nucleon and the pion cannot be at 

rest. Both have non-zero but opposite spatial momenta determined by the spatial volume 
and the boundary conditions imposed in the spatial directions. Still, for sufficiently large 
spatial volumes the exponent Mi can be smaller than the hrst one-particle excited state, 
associated with the Roper resonance A^*(1440) in infinite volume. Moreover, near physical 
quark masses the three-particle state containing the nucleon and two pions at rest will have 
a smaller energy than the one-particle excited state. 

The way to deal with multi-particle states in spectroscopy calculations is well-known. 
The well-established variational method [6] can be used provided interpolating fields for the 
multi-particle states are taken into account.^ Still, the more states one takes into account 
the larger is the generalized eigenvalue problem one has to solve numerically, and the error 
bounds for the energies derived in [10] get worse the denser the spectrum is. 

In this paper we provide some analytical results concerning the nucleon two-point func¬ 
tion. As has been pointed out in Ref. m, chiral perturbation theory (ChPT) can be 
employed to obtain an estimate for the ratio bi/bo-^ Moreover, to leading order in the chiral 
expansion one expects the ratio bi/bo to be independent of the a priori unknown low-energy 
constant (LEG) associated with the particular choice for the interpolating held. In that 
sense LO ChPT makes a rather dehnite prediction for bi/bo. Even if this result will receive 
substantial higher order corrections we do obtain a reliable hrst estimate for the impact of 
the nucleon-pion-state contribution to the two-point function. 


^ Very recent studies that include two-particle nucleon-pion states in the analysis are reported in Refs. [iH], 
for example. 

^ To our knowledge the idea for using ChPT to study the two-particle-state contributions to nucleon corre¬ 
lation functions was put forward first in Ref. |12j . 
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The results we find for 61/60 are small. For example, for a pion mass satisfying 4 

and ^ 0.2 we find 61/60 ~ 0.1. Taking into account the additional exponential 

suppression the two-particle-state contribution in ( 1 . 1 ) contributes at the few-percent level 
for euclidean times of about 0.5 fm. For larger t and in the effective mass the contribution 
is even smaller. Whether it is noticeable in practice is then a question of the size of the 
statistical errors in the lattice data. 

The nucleon-pion-state contribution to the two-point function has been independently 
computed in Ref. [I3]. The computation in that reference is performed in heavy baryon 
chiral perturbation theory (HBChPT), while here we employ the covariant formulation. In 
particular the chiral expressions for the interpolating fields differ in these two formulations, 
thus the final results are not the same. However, performing the appropriate expansion of 
our results we reproduce the results given in Ref. 


II. NUCLEON TWO-POINT CORRELATORS IN QCD 
A. General considerations 


Throughout this article we consider QCD in a finite spatial box. L denotes the box 
length in each direction and periodic boundary conditions are assumed. The euclidean time 
extent, however, is taken infinite. This choice implies an exponential decay of two-point 
functions and simplifies our calculations. Still, this simplification is a good approximation 
for many lattice QCD simulations. Another simplification concerns the masses of the up 
and down type quarks which we assume to be equal. Consequently, all three pions as well 
as the nucleons (proton and neutron) are mass degenerate. 

We are interested in the two-point correlation functions of a nucleon interpolating field 
N with definite parity. 


C±{t) = / d^x{N±{x,t)N±{0,0)) 


Il3 


(2.1) 


Here we defined N± = T±N and N± = NT± with the standard projectors F± = (70 ± l)/2.^ 
N itself is an interpolating field with the quantum numbers of the nucleon. Various choices 
are possible and we discuss concrete examples in the next subsection. For the moment we 
do not need to specify N any further. 

Let us consider the positive parity correlator. The integration over the spatial volume in 


Eq. (2.1) projects on states with zero total momentum. Hence, in the spectral decomposition 


of the correlator for large euclidean times t 3> 0 the dominant contribution comes from the 
single-particle state with the particle being the nucleon at rest. 


C+,N{t) = 


1 


2M4 


|(0|iV4(0)|iV(p = 0))|^e 


2 —Mn^ 


( 2 . 2 ) 


Here |V(p = 0)) is the nucleon state and Mjy denotes the nucleon mass. 

The interpolating field excites other states with the same quantum numbers as well. The 


contribution of an excited nucleon has the same form as Eq. (2.2) with the appropriate mass 


Note that in order to keep the notation simple we often suppress the Dirac and/or flavor indices. For 
example, eq. (2.1) contains an implicit summation over the Dirac indices. 
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M' > Mm- In addition we expect contributions from multi-hadron states. For sufficiently 
small pion masses the dominant multi-hadron states are those containing additional pions. 
For the two-particle nucleon-pion state contribution one finds 

■ (2.3) 

p ^ 


Etot is the total energy of the state and E,, are the individual energies of the nucleon 
and the pion, respectively. For weakly interacting pions F^tot equals approximately the 
sum Em + E^^. The sum over momenta runs over all momenta allowed by the boundary 
conditions imposed for the finite spatial volume, e.g. p = 27rn/L with n having integer¬ 
valued components. 

A simple dimensional analysis can be employed to make the volume suppression more 
quantitative. Assume the interpolating fields are local 3-quark-operators without derivatives. 


In that case the mass dimension of the matrix elements in (2.2) and (2.3) are 7/2 and 
5/2, respectively. Making the naive assumption (0|iV±(0)|A^(p)7r(—p)) (0|A^±(0)|iV±(p = 
0 ))//^ we can estimate the ratio of the two-particle and one-particle contributions as 


C+,V7r(t) ^ 1 

C'+,7v(t) ^ 2{ULfM^L E^ Em 


M„ Mm 


(2.4) 


If we assume the values ^ 200 MeV and L 4 fm we roughly find ^ 

1/30. The additional factors suppress the two-particle state contribution further, so we 
expect its contribution to the correlator to be rather small. 


States with more than one pion contribute analogously to (2.3), but each additional pion 
contributes an additional factor [ 2 (/ 7 rF)^MjrL]“^, i.e. the more pions in the state the larger 
the suppression of its contribution with the spatial volume. 

Our discussion applies to the negative parity correlator as well. In this channel the 
lightest single-particle state is, in infinite volume, the iV*(1535). However, the state with 
the lowest energy is the nucleon-pion state with both particles at rest, provided the pion 
mass is sufficiently small. Thus the two-particle state dominates the long time behaviour 


and the ratio analogous to the one in (2.4) will diverge for f —)■ cx) in the negative parity 
channel. 


B. Interpolating fields for the nucleon 

As already mentioned, there exist many choices for the interpolating field (“operator”) 
N with the quantum numbers of the nucleon. The number is significantly reduced if we 
consider local operators composed of three quark fields at the same point x. If, in addition, 
we constrain ourselves to operators without derivatives there exist only five different ones. 
As a consequence of Fierz identities only two are independent [HIIS] (see also [16]) and we 
focus on those. In order to write them down it is convenient to introduce the quark field 
doublet q as 

q = q^C'j5{ia2) - (2.5) 

Here q = (u, d)^ is the isospin doublet of the quark fields, C denotes the Dirac spinor charge 
conjugation matrix satisfying yj = —and a 2 is the second (isospin) Pauli matrix. 
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With these dehnitions the two nucleon operators can be written as 


Ni = {qq)q, 

N 2 = {.qi^qhbq- 


( 2 . 6 ) 


This compact form suppresses the contraction of the isospin and Dirac indices in the bilinear 
quark helds {qq) and {qj^q) (“diquarks”) and the summation over the color indices with an 
eabc to form a color singlet. The nucleon operators W are still isospin doublets. To project 
onto the quark content of the proton and neutron we need to contract with the isospin basis 
vectors Cp = (1, O)"*" and = (0,1)"*", respectively. However, in our case with preserved 
isospin symmetry any unit vector would be equally good. 

In the next section we need the counterparts of the nucleon operators W in ChPT. The 
mapping follows the standard procedure and rests on the transformation properties of W 
under chiral and parity transformations. The transformation properties under (singlet and 
non-singlet) chiral transformations have been studied in detail in Ref. [16]. Here we simply 
summarize the relevant results. 

We decompose the quark helds into right- and left-handed components, q = qR + qL, with 
the usual chiral projectors P+ = (l-|- 75)/2 and P_ = (1 — 75 ) 72 . It then follows that the held 
in (2.5) also decomposes according to q = qR + qi, with qR = qP+ and qR = qP-. The group 


of non-singlet chiral transformations is G = SU{2)r (g) SU{2)l, and under transformations 
L E G the chiral quark helds transform according to 


q = qR + qL RqR + LqR, 

i = in + it in -R* + 9l i'. 


(2.7) 


The diquarks decompose into qq = qRqR+qiqL and q%q = qRqR—qiqL, hence they transform 
as singlets under chiral transformations. Consequently, the transformation behavior of the 
nucleon operators is determined by the third quark held contribution q and 75 ^, given in 
(2.7). Decomposing this quark held into right- and left-handed components the complete 


nucleon helds Wg itself can be written as a sum of a right-handed and a left-handed term 
with the following transformation behavior under chiral transformations: 


Ni = N, 


R,L 


i,R 


+ Ni,L — RNi,R + LNi,L 


( 2 . 8 ) 


Concerning parity one hnds that qq and qj^q transform as a scalar and a pseudo scalar, 
respectively. Thus, both Ni and W transform as a Dirac spinor under parity, W —)■ 'joNi. 


Ni — Ni^R + Ni^R — > ''fo{Ni^R + Ni^r) 


(2.9) 


So far we considered local interpolating helds only. In lattice QCD so-called smeared 
interpolators are very often used, mainly to suppress excited-state contributions in the cor¬ 


relation function. Smeared nucleon interpolating helds are formed as in (2.6) but with the 
local quark helds replaced by smeared ones, which are generically of the form^ 


= / d‘}/K(x - y)q(y) 


( 2 . 10 ) 


^ We use a continuum notation here. In lattice QCD the integral is replaced by a sum over the lattice 
points. 
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with some gauge covariant kernel K{x — y) which is essentially zero for \x — y\ larger than 
some “smearing radius” R. The kernel depends on the details of the smearing procedure. 
Gaussian and exponential smearing dZHin] is local in time and the kernel contains a delta 
function in the euclidean time coordinate. In contrast, the gradient flow [20] is a truly 
four-dimensional smearing. 

What matters here are the transformation properties of the smeared quark helds. Pro¬ 
vided the kernel is diagonal in spinor space (as it is for Gaussian smearing and the gradient 
flow) the smeared quark helds transform just as the unsmeared ones under parity and global 
chiral transformations. Gonsequently, also the nucleon interpolating helds formed with the 
smeared quark helds transform according to (2.8) and (2.9), just as their local counterparts. 
Since the symmetry properties of the interpolating helds essentially determine their expres¬ 
sion in GhPT we can already conclude that both local and smeared interpolating helds are 
mapped onto the same ehective operator, dihering in their values for the LEGs only. We 
come back to this issue in section IIII Bi 

The two interpolating helds in eq. (2.6) were originally discussed in the context of QGD 
sum rule calculations, and in lattice QGD simulations the discretized version of these contin¬ 
uum interpolators are used. An alternative approach for the construction of baryonic opera¬ 
tors starts directly from the irreducible representations of the cubic group of the space-time 
lattice [2T|. In order to discuss this type of operators one hrst has to perform a mapping 
to the Symanzik ehective theory, the leading part being continuum QGD followed by cor¬ 
rections proportional to powers of the lattice spacing [22|. To the lattice operators whose 


leading Symanzik term is given by the interpolating helds in (2.6) the following discussion 
equally applies up to corrections proportional to the lattice spacing.® 

Finally, non-relativistic interpolators can be constructed from the relativistic ones if the 
standard non-relativistic (Dirac-)representation for the y-matrices is used 123 ]. In appendix 
[A] we discuss briehy the correlation function of the non-relativistic limit of Ni, which is 
rather easily obtained from the result using the relativistic Ni. 


III. THE NUCLEON TWO-POINT CORRELATORS IN CHPT 
A. The chiral Lagrangian 

The framework for our calculations is covariant baryon chiral perturbation theory pTl 
[23].® In this section we summarize a few relevant formulae since we work in euclidean space 
time and most references assume the Minkowski metric. 

We consider the chiral ehective Lagrangian^ 

= + (3.1) 

Here is the standard two-havor mesonic chiral Lagrangian to leading order mm- 


® Lattice artifacts are of course also present in lattice simulations that use the discretized expressions of Ni 
or Ay. 

® A thorough and pedagogical introduction to the subject can be found in Ref. (2^, for example. 

^ The superscripts denote the low-energy dimensions of these lagrangians, i.e. they count the number of 
derivatives and the power of quark mass terms |24j . 
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T u 

Ufj. 

R®L 

RT RuK^ = KuL^ 


P 

7oT ttl 



TABLE 1: Transformation behavior of the nucleon and pseudo scalar fields under chiral and parity 
transformations (see Ref. |29jl. The SU(2) matrix K appearing in the first row is defined by the 
transformation law of u such that = U transforms in the standard way. 


According to the conventions used here it reads 

4? = + t^TrlA1(l7 + C/t)!. ( 3 . 2 ) 

/, B are the standard LO LECs related to the pion decay constant and chiral condensate in 
the chiral limit.® The pion helds are contained in the held U according to 


U (x) = exp 


yvr [x}a 


(3.3) 


with the usual Pauli matrices cr“. Ai denotes the quark mass matrix. With equal quark 
masses m for the up and down quark it is proportional to the unit matrix. In that case all 
three pions have the same mass which to LO is related to the quark mass via = 2Bm. 


The second part in the chiral lagrangian (3.1) contains the nucleon helds and their cou¬ 
pling to the pions, 

= ^{U) + Mn- (3.4) 

The helds T = (p, n)^ and T = (p, n) denote the nucleon helds with two Dirac spinors for 
the proton p and the neutron n. and qa are the nucleon mass and the axial-vector 

coupling constant in the chiral limit. Since we assume isospin symmetry the proton and the 
neutron are mass degenerate. 


The pion helds enter via the held the so-called chiral vielbein, dehned by 


u, 


= i[u^d^u — ud^u^], u{x) = \/U{x). 


(3.5) 


A second source of pion-nucleon coupling stems from the covariant derivative ij) = 


m 


(3.4), with -b T^^T and T^ = [u^di^u + udf,u^ /2. 

The construction of the chiral lagrangian is based on the symmetry properties of the 
underlying QCD lagrangian, which Ces needs to reproduce. The transformation behavior 
of the nucleon and pseudo scalar helds under chiral and parity transformations is briehy 
summarized in table (for details see Ref. [26], for example). 


Expanding and T^ we obtain pion-nucleon interaction terms with various numbers of 
pion helds. Since is parity-odd and T^ is parity-even the leading interaction term with 
one pion held only stems from Un and reads 


4.1lo = 


(3.6) 


® Our conventions correspond to /,r = 92.2 MeV. 
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This interaction term couples two axial vectors to obtain a Lorentz scalar. In addition, 
isospin symmetry is preserved. 

The chiral Lagrangian incorporates a derivative expansion and the chiral dimension counts 
the number of derivatives and powers of the quark mass. The complete list of terms through 
fourth order is given in [30] • For the purpose of this paper, however, the term in (3.6) is 
sufficient. 


For the perturbative calculation in section [III C| we need the propagators for the nucleon 
and the pion in position space. The pion propagator is the same as in Ref. El. 




1 


^ipix-y) ^-E^lxo-yol 


(3.7) 


with pion energy + M^. The nucleon propagator S'“^(x,|/) is also easily derived 

from the quadratic term in (3.4), and it reads 




= 5 “"^ 


ab T —3 


2 ’^ 

^ip(x-y)^-E,^\xo-yo\ 


2E 


N 


(3.8) 


a, b and af3 refer to the isospin and Dirac indices, respectively. The factor Zp (spinor indices 
suppressed) in the numerator is defined as 

Zp = -ip- 7 ± EnJo + Mn , (3.9) 

where the + (—) sign applies to xq > yo (xq < yo), and the nucleon energy En = + Mfj 

involves the nucleon mass. The sum in both propagators runs over the discrete spatial 
momenta that are compatible with periodic boundary conditions, i.e. p = 27rf?/L with n 
having integer-valued components. 


B. The chiral expansion of the interpolating fields 


The construction of the nucleon operators in baryon ChPT follows the standard proce¬ 
dure. Based on the symmetry properties of the operators on the quark level we write down 
the most general expression in the effective theory that has the same symmetry properties. 
This has essentially been done in Ref. [29] and we summarize the results needed in the 
following.® 


The nucleon operators in the effective theory needs to transform as given in (2.8) and 


(2.9) under chiral and parity transformations. Basically, the nucleon operators are a sum 
of a right- and left-handed spinor and can be written as (we follow closely the notation 
introduced in Ref. 


N = 


^ ^ 




a 


At'- -I- N 


(n) 


k,L 


(3.10) 


® Instead of IVi _2 two other operators, related to IVi 2 by Fierz identities, are considered in Ref. jH] . However, 
the chiral expansion is essentially the same. 











n 

k 

^^k,R 


0 

1 



1 

1 

uUfj,P+'y^A> 


1 

2 

uu^P+D^Af 

-u^Uf,P-D^Ai 


TABLE 2: Low-energy operators for the nucleon interpolating fields through chiral dimension one 
(see Ref. [29]). 


and are operators with low-energy dimension n. in denotes the nnmber of oper- 


'fc,L 


ators with chiral dimension n, which are labelled by the index k. Under chiral and parity 


transformations the fields in (3.10) transform according to 


r(’^) 


Nn + N, 


(n) 

k,L 

(n) 


A'tX + Kl 


P / I ArU)'i 
;,L ^^k,R) ■ 


-)■ 


7o(A'r, 


( 3 . 11 ) 


Each term on the right hand side of (3.10) comes with its own LEG and it is parity 
that relates the coefficients of the right- and left-handed contribntions. 

An incomplete list of operators throngh chiral dimension two can be fonnd in Ref. [2^ . 
For convenience we reproduce the ones through n = 1 in table It is straightforward to 
check that these operators satisfy the transformation laws given in (3.11). 


Note that N in (3.10) does not carry an index i that would refer explicitly to one of 
the two interpolating fields defined in (2.6). We dropped this index because the chiral 
expansion for both operators is the same due to their similar transformation behavior. The 
only difference are different values for the LECs in the chiral expansion. In order to keep 
our notation simple we suppress an additional index at the operator and the LECs in the 
following. 


Similarly, for the smeared interpolating fields discussed at the end of section IIB we also 


find the same effective operator (3.10) with different LECs. However, one qualification has 
to be made. Smeared interpolators with some “size” are mapped onto the pointlike nucleon 
field in the chiral effective theory. For this to be a good approximation the smearing radius 
needs to small compared to the Compton wave length of the pion. Provided this condition 
is met the pions do not distinguish between smeared and pointlike interpolating fields. 

At lowest low-energy dimension only one operator contributes. Expanding u, in powers 
of pion fields and keeping only the terms up to linear order we obtain for N the expression 


N{x) = a ( T(a:) -F —7r{x)'y5'^{x) 


a = Aa 


( 0 ) 


( 3 . 12 ) 


The first LO term is proportional to the nucleon field T, as expected. The second NLO term 
(suppressed by 1//) involves a nucleon-pion coupling that will contribute to the two-particle 
nucleon-pion terms (|2.3) in the nucleon correlation function. 
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A concrete example is given in Ref. m where ChPT for some observables based on the gradient flow has 


been constructed. 


9 



















£1 


a) 



b) 



c) 



• 4 


FIG. 1: Feynman diagrams for the nucleon correlation function. The squares represent the nucleon 
operator at times t and 0, where the open and solid squares denote the leading and next-to-leading 


order terms given in (3.12). The circles represent a vertex insertion at an intermediate space time 


point; and an integration over this point is implicitly assumed. The solid and dashed lines represent 
nucleon and pion propagators, respectively. 


C. Perturbative expansion of the correlation functions 


We are now in the position to compute the correlation functions (2.1) perturbatively 
within the chira l effec tive theory. The leading contribution is obtained by taking into account 
the LO term in (3.12) for N and N = iV^yo. Since these helds are proportional to the nucleon 
helds T, T the LO contribution is essentially the nucleon propagator. In terms of Feynman 
rules in position space this contribution is represented by the Feynman diagram in hgure 
[J. Taking into account ( |3.8[ ) for the nucleon propagator the LO results for the correlators 
are easily obtained, 


C+,N{t) = 2|d|V^^*, = 0. (3.13) 

These results are a single-particle state contributions to the correlation function, and by 
comparing with (2.2) we can read off the LO relation between the vacuum-to-nucleon matrix 
element and the LEG a, 


|{0|]V+(0)|]V+(p = 0))|" = 4M„|a| 


(3.14) 


in (3.13) vanishes at this order because our effective theory does not contain the 
negative parity nucleon as a degree of freedom. 

The diagrams in Figs. m-in form the leading contribution to the two-particle nucleon- 
pion part of the correlation function. The dashed line represents the pion propagator, which, 
together with the nucleon propagator, leads to terms with the expected exponential fall-off 
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n 0 

1 2 

3 

4 5 6 

7 8 

rUn 1 

6 12 

8 

6 24 24 

0 12 


TABLE 3: Multiplicities in eq. (3.17) for n < 8 (see Ref. [32]). 
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with Etot = En + E^ 

The calculation of the diagrams [^ - [^ is straightforward, and the final results can be 
compactly written as 


C+^Nnit) = 2 |Q!|^ 
C-^Nnit) = 2 | q ;|^ 


3 V ^ TTl E]\j — -^N 
8 (/L) 2 mL^^ 


E 


2 i?jv 
m En + Mat 


8 (/L) 2 mL ^ E^ 

p 


2E 


N 


9a 


9a 


Etot + Mj 


N 


Etot ~ M, 


N 


-,-Eu 


, (3.15) 


Etot ~ M, 


1 2 


N 


Etot + 


N 


3— -E^tot 1^1 


(3.16) 


Note that the p = 0 term in C+^ 7 V 7 r(^) correctly vanishes, as required by parity. 

The overall factor 2\a\^ has its origin in the appearance of lap as an overall factor for the 
two terms in (3.12). This implies that the relative size of the two-particle state contributions 


in (3.15), (3.16) and the one-particle state contribution in (3.13) does not contain the LEG 


d associated with the nucleon interpolating field. It only depends on the LECs / and qa of 
the effective action. However, this is true to LO only. Taking into account the higher order 
operators in table their LECs will not cancel in the ratio. 

Recall that the chiral expansion for the nucleon operators Ni and N 2 are the same, the 
only difference being different LECs. Restoring the label i = 1, 2 in the LEG d the results 
given above (with |dp replaced by |djp) refer to the correlation functions where the same 
operator is used for both source and sink. In case of the correlator with Ni{x) and iV2(0) 
we simply need to replace |dp by did^. Since the same combination appears as an overall 
factor in the single-particle and two-particle state contributions it still drops out in the ratio 
of these two contributions. The same statement applies to smeared operators. 


The summation in (3.15), (3.16) is over all lattice momenta compatible with periodic 
boundary conditions. Those that are related by the symmetries of the spatial lattice lead 
to the same contribution, hence it is convenient to sum over the absolute value p = |p|. 
Imposing periodic boundary conditions the absolute value can assume the values p„ = 
[2tt/L) t/n, n = n\ + n\ + with the nk being integers, and the sums in the results given 
above are replaced according to 

(3.17) 


E^E 


rrir. 


Pn 


The multiplicities count the number of vectors p with the same Pn- Multiplicities for 
n < 20 are given in Ref. [32] (for convenience we summarize the hrst eight in table |^. 


Fig. & also contains a contribution with a time dependence proportional to texp(—MArt). This results in 


the renormalization of the nucleon mass and can be ignored for our purposes. 
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D. Final Results 


Adding the two results in (3.13) and (3.15) the positive parity correlation function can 
be written as 


C+{t) = 2|d| 




g~ (-^tot, n —t 


Pn 


where we introduced new dimensionless coefficients 

L ^rur,. 


8{fL)^E^,nL 


K, 


ht = 


EN,n — Af, 


N 


2E 


N,n 


9a 




tot,n 


+ M^ 


N 


Etot,n ~ Af, 


N 


(3.18) 

(3.19) 

(3.20) 


Result (3.18) corresponds to the example we discussed briefly in the introduction. The 
coefficients are equal to bn/bo and all nucleon-pion state contributions are given. 

We already mentioned that the positive parity two-point correlator was independently 
calculated in Ref. PI using HBChPT |33l |3l]. In that non-relativistic formulation one 
drops the antinucleon degrees of freedom and the dispersion relation of the heavy nucleon is 
non-relativistic. If we expand 

terms we hnd g\p'^/E^ and reproduce the result in 

The negative parity channel is slightly different since there is no single particle state 
contribution stemming from the nucleon. The leading single-particle state contribution 
comes from the negative parity partner N* which is not a degree of freedom in our effective 
theory. And even if we included it explicitly the coupling of the interpolating field N~ to 


Mn+pI^/2Mn in (3.20) and drop all but the dominant 


the N* would come with a LEG unrelated to the LEG a that enters result (3.16). 


The dominant contribution in our result for C- stems from the nucleon-pion state with 
the nucleon and the pion at rest. Taking this contribution out of the sum we arrive at the 
form 


C-{t) 


2|^|2 _ p-iMN+MT,)t 

' ' 8{fL)m^L 


^ ^ Q g (-^tot,n -A^7r)i 


(3.21) 


with a coefficient 


= 


h„ = 


A- 

m „,—— K 


Z? 

^77,n 

EN,n + Afj 


N 


2E 


N,n 


9a 


Etot,n — M, 


N 


Etot,n + Af, 


N 


(3.22) 

(3.23) 


Note that the results for the two parity channels are proportional to the unknown LEG 
I dp. Thus, taking the ratio G_(f)/G+(f) this constant drops out and the only LEGs con¬ 
tributing to this ratio are / and qa of the LO chiral lagrangian. 


E. Numerical estimates 

Before trying to estimate the impact of the nucleon-pion-state contribution to the two- 
point function it is necessary to discuss the conditions for the applicability of the results 
derived in the last section. 
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A4 

M^L 

M^/Ay Pn/Ay 

Umax 

200 

4 

~ 1/6 ~ yJn/A: 

1 

150 

4 

~ 1/8 ~ \/n/5 

2 

130 

4 

~ 1/9 ~ \/n/6 

3 

150 

6 

~ 1/8 ~ \/n/8 

5 


TABLE 4: The chiral expansion parameters and pn /various pion masses and spatial 

extensions L. Umax in the last column stems from the condition Pn^ax/'^x ~ 0-3 


ChPT is an expansion in the pion mass and momentnm. Both need to be small compared 
to the chiral symmetry breaking scale which is typically identihed with 47r/^. In a finite 
spatial volnme with periodic bonndary conditions the pion momenta are discrete and the 
condition for the applicability of the chiral expansion reads [32] 


A. 




n 


2U M^L 


< 1 . 


(3.24) 


Even though the pion mass cancels on the right hand side we prefer this form since lattice 
QCD configurations are often characterized in terms of the pion mass and M^^L. Given these 
two numbers and the pion decay constant ~ 90MeV eq. (3.24) provides a bound on the 


pion momenta and the label n. Table lists a few representative values that approximately 
match the parameters in present-day lattice simulations.^^ The expansion parameter M^/A^ 
is sufficiently small for pion masses smaller than 200MeV. The situation is less favourable 
for pn/Ay.. In order to have at least theoretically a chiral expansion at all the expansion 
parameter Pn/Ay should be reasonably smaller than 1. If we restrict ourselves to momenta 
satisfying Pn/Ay < 0.3 we find values for Umax ranging between 2 and 5. For the coefficients 

with n < Umax we expect the chiral expansion to be applicable. 

Besides the question of applicability there is the question of how rapidly the chiral ex¬ 
pansion converges. The smaller n and Pn is the better the chiral expansion behaves. For 
Pn!Ay ~ 0.3 one does not expect a fast rate of convergence, and the LO results for the coeffi¬ 
cients associated with these rather high momenta will probably receive rather large higher 
order corrections. Definite statements about these higher order corrections are difficult to 
make without having done the calculation, but as a rough error estimate we may allow for 
a 50% error. The error will be smaller for the coefficients associated with the smaller mo¬ 
menta. Note that the exponential suppression due to the exponential exp[—(Ftot.n — AI]x)t] 
is stronger for the contributions with larger momentum, so the contributions with a larger 
uncertainty in the two-point function are more suppressed. 

Another reason to constrain the momenta Pn and n stems from the requirement that the 
nucleon-pion-state energy should be sufficiently well separated form the first one-particle 
resonance energy. If that is not the case one expects large mixing between these two states 


For the calculation of the light hadron spectrum in [3S] the BMW collaboration generated a lattice 
ensemble with « 190MeV with « 3.9. The PACS collaboration has recently reported results 
obtained on a 96“^ lattice with ~ 147MeV and « 6 [36]. Finally, the FERMILAB and MILC 
collaborations [3] have generated a lattice ensemble with Mt^ « 128MeV and Mt^L « 3.9. 
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FIG. 2: as a function of the sink-source separation t for the parameter sets in 


table 


that significantly alter the coefficents cjJ 


In practice, however, this seems to be almost 


automatically satished once the momenta are constrained by the bound (3.24), because the 
energy of the hrst resonance is about O.SGeV higher than the nucleon mass. For example, 
for the values Umax given in table the total energy of the nucleon-pion state is at least 
lOOMeV smaller than the expected resonance energy, and the energy gap for the states with 
n < n-max is even smaller. We therefore expect mixing effects to be negligible, at least within 
the uncertainties the LO results presented here are afflicted with. 

In lattice simulations one usually computes the effective nucleon mass, dehned as the 
negative time derivative of logG+(f). With (3.18) we obtain 


MN,eS = M, 


dt = 


N 


1 




Pn 




tot,n 


M, 


N 


(3.25) 

(3.26) 


Figure]^ shows the ratio Mn^ss/Mn as a function of t for the four parameter sets in table]^ 
The deviation from the constant value 1 is caused by nucleon-pion-state contribution. The 
sum over pn is truncated at Umax given in table We conclude that the nucleon-pion-state 
contribution to the effective mass is small, less than 2 percent for t > 0.5fm, and dropping 
below 1 percent at f ~ Ifm. Note that the nucleon-pion-state contribution to the effective 
mass is signihcantly smaller than to the two-point function itself since the coefficients are 
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smaller than the coefficients c^. For the parameter sets of table the factor Etot,n/MN — 1 
cansing this snppression is less than abont 0.4. 

Fignre shows the resnlts to leading order in the chiral expansion. As argned before 
we may allow for a 50% error for the coefficents, which approximately resnlts in a 50% 
error in hgure Even with this nncertainty we can conclnde that the nncleon-pion-state 
contribntion to the effective mass is at the few-percent level. Whether this plays a role in 
practice depends on the statistical errors in the lattice data and the sonrce-sink separations 
accessible in the simnlation. For statistical errors of 1% and below the nncleon-pion-state 
contribntion shonld be noticeable in the data, at least for enclidean times below 0.7 fm. 
Obvionsly, actnal lattice data needs to be analyzed with the resnlts presented here before 
dehnite conclnsions can be drawn. 


IV. CONCLUDING REMARKS 

The results we hnd for the nncleon-pion-state contribution to the nucleon correlation 
function is small. Whether it is too small to be seen in lattice QCD data is not easy to 
answer and depends on various parameters of the simulation, among others on the range for 
the euclidean time and the size of the statistical errors in the data. Future physical-point 
simulations aiming at a one-percent error for the nucleon mass will probably be sensitive to 
the corrections studied here. 

We repeatedly mentioned that the calculations in this paper are LO calculations. In 
principle, the higher order corrections can be calculated straightforwardly. Whether this is 
useful in practice is somewhat doubtful. The NLO contributions will depend on the LECs 
associated with the interpolating nucleon helds. Since their values are a priori unknown the 
NLO results are in a way less predictive than the LO results derived here. On the other 
hand, the LO results are universal for all the various nucleon interpolating helds considererd 
here. In order to discuss and describe differences caused by using different interpolating 
helds one has to go at least one order higher in the chiral expansion. 

Here we studied the dominant nncleon-pion-state contributions to the nucleon two-point 
function. The same framework can be used to study the nncleon-pion-state contribution to 
higher n-point functions. The nucleon axial charge, for example, is obtained from the three- 
point function involving the axial vector current. In that case the contamination with excited 
and multi-particle states is expected to be larger, because more than one euclidean time 
diherence are present in these correlation functions, and these are usually smaller than the 
source-sink separation in the two-point function. Preliminary results concerning the nncleon- 
pion-state contributions in the determination of qa can be found in [32] • Other interesting 
observables to study are the electromagnetic form factors and the quark momentum fraction 
in the nucleon, for instance. 
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We mention that Ref. [38] reports on the presence of a non-negligible nucleon-pion-state contribution to 
C-. Unfortunately, a direct check of our results is not possible. 
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Appendix A: Non-relativistic interpolating fields 


As mentioned at the end of section IIB non-relativistic interpolating nucleon fields can 
be constructed from their relativistic counterparts. Here we briefly discuss the minor mod¬ 
ification this choice implies in case of Ni. It is sufficient to discuss the local interpolating 
field only, the generalization to the smeared fields is completely analogous to the discussion 
at the end of section Ell 

In our discussion of the relativistic interpolating fields there was no need to specify a 
particular representation for the 7-matrices. For the non-relativistic limit it is useful to 
assume the Dirac representation with 70 being diagonal. The reason is that in this basis the 
projectors F^ = (70 ± l)/2 can be used to decompose a quark spinor q into its upper (large) 
spinor component and its lower (small) spinor components qi. The latter vanish in the 
non-relativistic limit. 

Decomposing Ni into upper and lower components one finds that one contribution in¬ 
volves upper components only ra. Let us call this part Ai^nr = (?«?«)?«• This field com¬ 
pletely breaks charge conjugation but it still transforms according to (2.8) and (2.9) under 
chiral transformations and parity. So this field is mapped onto the expression (3.12) in 
the chiral effective theory but the nucleon interpolating field T in the terms of table |2 are 
replaced by F+T. Using this we find the expression 




[X = Or, 


F+T(a:) - T.—7i{x)^5'^{x) 


(Al) 


for the interpolating field in the chiral effective theory. In addition to a different LEG the 


parity projectors appear explicitly in this expression, in contrast to the result in (3.12). 


Suppose we compute the positive parity correlation function with this interpolating field. 
Due to the parity projector F_|_ in the definition of the correlator only the first term in (Al) 


contributes, the second vanishes identically. In terms of Feynman diagrams this means that 
only diagram e) in fig. contributes. Dropping the contributions from the other diagrams 


in the result (3.15) is very simple, it amounts in dropping the 1 in [1 — g^.. .]^ on the right 


hand side of (3.15). This is esily understood since diagram e) is the only one with two vertex 
insertions and therefore the only one proportional to g\. This is the only modification there 
is, thus the final result for the correlator is the same as in (3.18) provided the aforementioned 
modification is made in the coefficient h^. Note that dropping the 1 in [1 — .. .]^ makes 

the coefficient larger, even though the increase is not huge, at least for small energies 
where the contribution proportional to gA is much larger than 1. 

A similarly simple modification is found in the negative parity correlator. In that case 


only the second term in (Al) provides a non-vanishing contribution to the correlator, thus 


only diagram b) contributes. This is accounted for by keeping only the 1 in [1 — ...] 


2 ; 


m 


( |3T6| and ( |3^ , respectively. 
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